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^ ! Abstract 
<D 

A measurement of the forward-backward asymmetry of e + e~ — > cc and e + e~ — > 
bb on the Z resonance is performed using about 3.5 million hadronic Z decays 
collected by the DELPHI detector at LEP in the years 1992 to 1995. The heavy 
quark is tagged by the exclusive reconstruction of several D meson decay modes. 
The forward-backward asymmetries for c and b quarks at the Z resonance are 
^ ■ determined to be: 

A c FB {s/s = 91.235 GeV) = 0.0659 ± 0.0094 (stat) ± 0.0035 (syst) 

A h FB {</s = 91.235 GeV) = 0.0762 ± 0.0194 (stat) ± 0.0085 (syst) 

A c fb(\/s = 89.434 GeV) = -0.0496 ± 0.0368 (stat) ± 0.0053 (syst) 

Afb(Vs = 89.434 GeV) = 0.0567 ± 0.0756 (stat) ± 0.0117 (syst) 

A c FB (^s~ = 92.990 GeV) = 0.1180 ± 0.0318 (stat) ± 0.0062 (syst) 

Afb(V~s = 92.990 GeV) = 0.0882 ± 0.0633 (stat) ± 0.0122 (syst) 

The combination of these results leads to an effective electroweak mixing angle 
of: 

sin 2 0gj = 0.2332 ± 0.0016 



(Submitted to Eur. Phys. J. C) 



ii 



P.Abreu 21 , T.Adye 36 , P.Adzic 11 , Z.Albrecht 17 , T.Alderweireld 2 , G.D.Alekseev 16 , R.Alemany 49 , T.Allmendinger 17 
P.P.Allport 22 , S.Almehed 24 , U.Amaldi 9 , S.Amato 47 , E.G.Anassontzis 3 , P.Andersson 44 , A.Andreazza 9 , S.Andringa 21 
P.Antilogus 25 , W-D.Apel 17 , Y.Arnoud 9 , B.A sman 44 , J-E.Augustin 25 , A.Augustinus 9 , P.Baillon 9 , P.Bambadc 19 

F. Barao 21 , G.Barbiellini 46 , R.Barbier 25 , D.Y.Bardin 16 , G.Barker 17 , A.Baroncelli 38 , M.Battaglia 15 , M.Baubillier 23 
K-H. Becks 52 , M.Begalli 6 , P.Bcillicrc 8 , Yu.Bclokopytov 9 ' 53 , K.Belous 42 , A.C.Bcnvcnuti 5 , C.Bcrat 14 , M.Berggrcn 25 

D. Bertini 25 , D.Bertrand 2 , M.Besancon 39 , F.Bianchi 45 , M.Bigi 45 , M.S.Bilenky 16 , M-A.Bizouard 19 , D. Block 10 
H.M.Blom 30 , M.Boncsini 27 , W.Bonivcnto 27 , M.Boonekamp 39 , P.S.L.Bootb 22 , A.W.Borgland 4 , G.Borisov 19 , C.Bosio 41 
O.Botncr 48 , E.Boudinov 30 , B. Bouquet 19 , C.Bourdarios 19 , T.J.V.Bowcock 22 , I.Boyko 16 , I.Bozovic 11 , M.Bozzo 13 
P.Branchini 38 , T.Brenke 52 , RA. Brenner 48 , P.Bruckman 18 , J-M.Brunet 8 , L.Bugge 32 , T.Buran 32 , T.Burgsmueller 52 
P.Buschmann 52 , S.Cabrera 49 , M.Caccia 27 , M.Calvi 27 , T.Camporesi 9 , V.Canale 37 , F.Carena 9 , L.Carroll 22 , C.Caso 13 
M.V.Castillo Gimcncz 49 , A.Cattai 9 , F.R.Cavallo 5 , V.Cbabaud 9 , M.Chapkin 42 , Ph.Cbarpentier 9 , L.Cbaussard 25 
P.Cbeccbia 35 , G.A.Cbelkov 16 , R.Chicrici 45 , P.Cbliapnikov 42 , P.Chocbula 7 , V.Chorowicz 25 , J.Cbudoba 29 , K.Cieslik 18 
P.Collins 9 , R.Contri 13 , E.Cortina 49 , G.Cosme 19 , F.Cossutti 9 , J-H.Cowell 22 , H.B.Crawley 1 , D.Crcnnell 36 , S. Crepe 14 

G. Crosetti 13 , J. Cucvas Maestro 33 , S.Czellar 15 , M.Davenport 9 , W.Da Silva 23 , A.Degborain 2 , G. Delia Ricca 46 
P.Delpierre 26 , N.Dcmaria 9 , A.Dc Angelis 9 , W.De Boer 17 , S.Dc Brabandere 2 , C.De Clercq 2 , B.De Lotto 46 , A.De Min 35 
L.Dc Paula 47 , H.Dijkstra 9 , L.Di Ciaccio 37 ' 9 , J.Dolbeau 8 , K.Doroba 51 , M.Dracos 10 , J. Drees 52 , M.Dris 31 , A.Duperrin 25 
J-D.Durand 9 , G.Eigen 4 , T.Ekelof 48 , G.Ekspong 44 , M.Ellcrt 48 , M.Elsing 9 , J-P.Engel 10 , B.Erzen 43 , M.Espirito Santo 21 

E. Falk 24 , G.Fanourakis 11 , D.Fassouliotis 11 , J.Fayot 23 , M.Feindt 17 , A.Fenyuk 42 , P.Ferrari 27 , A.Fcrrcr 49 

E. Ferrer-Ribas 19 , S.Fichet 23 , A. Firestone 1 , U.Flagmeyer 52 , H.Foetb 9 , E.Fokitis 31 , F.Fontanelli 13 , B.Franek 36 

A. G.Frodesen 4 , R.Fruhwirtb 50 , F.Fulda-Qucnzer 19 , J.Fuster 49 , A.Galloni 22 , D.Gamba 45 , S.Gamblin 19 , M.Gandclman 47 
C.Garcia 49 , C. Caspar 9 , M. Caspar 47 , U.Gasparini 35 , Ph.Gavillet 9 , E.N.Gazis 31 , D.Gele 10 , N.Ghodbane 25 , I.Gil 49 

F. Glege 52 , R.Gokieli 9 ' 51 , B.Golob 43 , G.Gomcz-Ceballos 40 , P.Goncalves 21 , I.Gonzalcz Caballcro 40 , G.Gopal 36 
L.Gorn 1 ' 54 , M.Gorski 51 , Yu.Gouz 42 , V.Gracco 13 , J.Grahl 1 , E.Graziani 38 , C.Green 22 , H-J. Grimm 17 , PGris 39 

G. Grosdidier 19 , K.Grzelak 51 , M.Guntber 48 , J.Guy 36 , F.Hahn 9 , S.Habn 52 , S.Haider 9 , A.Hallgren 48 , K.Hamacher 52 
J.Hansen 32 , F.J.Harris 34 , V.Hedberg 24 , S.Heising 17 , J.J.Hernandez 49 , P.Herquet 2 , H.Herr 9 , T.L.Hessing 34 
J.-M.Heuser 52 , E.Higon 49 , S-O. Holmgren 44 , P.J.Holt 34 , S.Hoorelbeke 2 , M.Houlden 22 , J.Hrubec 50 , K.Huet 2 
G.J.Hughes 22 , K.Hultqvist 44 , J.N.Jackson 22 , R. Jacobsson 9 , P.Jalocha 9 , R.Janik 7 , Ch.Jarlskog 24 , G.Jarlskog 24 
P.Jarry 39 , B.Jean-Marie 19 , E.K.Johansson 44 , P.Jonsson 25 , C.Joram 9 , P.Juillot 10 , F.Kapusta 23 , K.Karafasoulis 11 
S.Katsanevas 25 , E.C.Katsoufis 31 , R.Kcrancn 17 , B.P.Kersevan 43 , B.A.Khomenko 16 , N.N.Khovanski 16 , A.Kiiskinen 15 

B. King 22 , A.Kinvig 22 , N.J.Kjaer 30 , O.Klapp 52 , H.Klein 9 , P.Kluit 30 , P.Kokkinias 11 , M.Koratzinos 9 , V.Kostioukhinc 42 

C. Kourkoumclis 3 , O.Kouznctsov 39 , M.Krammer 50 , E.Kriznic 43 , P.Krstic 11 , Z.Krumstcin 16 , P.Kubincc 7 , J.Kurowska 51 
K.Kurvinen 15 , J.W.Lamsa 1 , D.W.Lane 1 , P.Langefeld 52 , V.Lapin 42 , J-P.Laugier 39 , R.Lauhakangas 15 , G.Leder 50 
F.Ledroit 14 , V.Lefebure 2 , L.Leinonen 44 , A.Leisos 11 , R.Leitner 29 , J.Lemonne 2 , G.Lenzen 52 , V.Lepeltier 19 , T.Lesiak 18 
M.Lethuillier 39 , J.Libby 34 , D.Liko 9 , A.Lipniacka 44 , I.Lippi 35 , B.Loerstad 24 , J.G.Loken 34 , J.H.Lopes 47 , J.M.Lopcz 40 
R.Lopez-Fernandez 14 , D.Loukas 11 , P.Lutz 39 , L.Lyons 34 , J.MacNaughton 50 , J.R.Mahon 6 , A.Maio 21 , A.Malek 52 
T.G.M.Malmgren 44 , V.Malychev 16 , F.Mandl 50 , J.Marco 40 , R.Marco 40 , B.Marcchal 47 , M.Margoni 35 , J-C. Marin 9 
C.Mariotti 9 , A.Markou 11 , C.Martinez-Rivero 19 , F. Martinez- Vidal 49 , S. Marti i Garcia 9 , J.Masik 12 
N.Mastroyiannopoulos 11 , F.Matorras 40 , C.Matteuzzi 27 , G.Matthiae 37 , F.Mazzucato 35 , M.Mazzucato 35 , M.Mc Cubbin 22 
R.McKay 1 , R.Mc Nulty 22 , G.Mc Phcrson 22 , C.Mcroni 27 , W.T.Meyer 1 , E.Migliore 45 , L.Mirabito 25 , W.A.Mitaroff 50 
U.Mjoernmark 24 , T.Moa 44 , M.Moch 17 , R.Moeller 28 , K.Moenig 9 , M.R.Monge 13 , X.Moreau 23 , P.Morettini 13 , G.Morton 34 
U.Mueller 52 , K.Muenich 52 , M.Mulders 30 , C.Mulet-Marquis 14 , R.Murcsan 24 , W.J.Murray 36 , B.Muryn 14 ' 18 , G.Myatt 34 
T.Myklebust 32 , F.Naraghi 14 , F.L.Navarria 5 , S.Navas 49 , K.Nawrocki 51 , P.Negri 27 , S.Nemecek 12 , N.Neufeld 9 
N.Neumeister 50 , R.Nicolaidou 39 , B.S.Nielsen 28 , M.Nikolenko 10 ' 16 , V.Nomokonov 15 , A.Normand 22 , A.Nygren 24 
V.Obraztsov 42 , A.G.Olshcvski 16 , A.Onofre 21 , R.Orava 15 , G.Orazi 10 , K.Osterberg 15 , A.Ouraou 39 , M.Paganoni 27 
S.Paiano 5 , R.Pain 23 , R.Paiva 21 , J.Palacios 34 , H.Palka 18 , Th.D.Papadopoulou 31 , K.Papageorgiou 11 , L.Papc 9 , C.Parkes 9 
F.Parodi 13 , U.Parzefall 22 , A.Passeri 38 , O.Passon 52 , M.Pegoraro 35 , L.Peralta 21 , M.Pernicka 50 , A.Perrotta 5 , C.Petridou 46 
A.Petrolini 13 , H.T.Phillips 36 , F.Pierre 39 , M.Pimcnta 21 , E.Piotto 27 , T.Podobnik 43 , M.E.Pol 6 , G.Polok 18 , P.Poropat 46 
V.Pozdniakov 16 , P.Privitera 37 , N.Pukhaeva 16 , A.Pullia 27 , D.Radojicic 34 , S.Ragazzi 27 , H.Rahmani 31 , D.Rakoczy 50 
P.N.Ratoff 20 , A.L.Read 32 , P.Rebecchi 9 , N.G.Redaelli 27 , M.Regler 50 , D.Reid 30 , R.Rcinhardt 52 , P.B.Rcnton 34 
L.K.Resvanis 3 , F.Richard 19 , J.Ridky 12 , G.Rinaudo 45 , O.Rohne 32 , A.Romcro 45 , P.Ronchese 35 , E.I.Rosenberg 1 
P.Rosinsky 7 , P.Roudeau 19 , T.Rovelli 5 , Ch.Royon 39 , V.Ruhlmann-Klcidcr 39 , A.Ruiz 40 , H.Saarikko 15 , Y.Sacquin 39 

A. Sadovsky 16 , G.Sajot 14 , J. Salt 49 , D.Sampsonidis 11 , M.Sannino 13 , H. Schneider 17 , Ph.Schwcmling 23 , U.Schwickerath 17 
M.A.E.Schyns 52 , F.Scuri 46 , P.Seagcr 20 , Y.Scdykh 16 , A.M.Segar 34 , R.Sekulin 36 , R.C.Shellard 6 , A.Sheridan 22 , M.Siebel 52 
L.Simard 39 , F.Simonctto 35 , A.N.Sisakian 16 , G.Smadja 25 , N.Smirnov 42 , O.Smirnova 24 , G.R.Smith 36 , A.Sopczak 17 
R.Sosnowski 51 , T.Spassov 21 , E.Spiriti 38 , P.Sponholz 52 , S.Squarcia 13 , D.Stampfer 50 , C.Stancscu 38 , S.Stanic 43 
K.Stevenson 34 , A.Stocchi 19 , J.Strauss 50 , R.Strub 10 , B.Stugu 4 , M.Szczckowski 51 , M.Szeptycka 51 , T.Tabarelli 27 
F.Tegenfeldt 48 , F.Terranova 27 , J.Thomas 34 , J.Timmcrmans 30 , N.Tinti 5 , L.G.Tkatchcv 16 , S.Todorova 10 , A.Tomaradze 2 

B. Tome 21 , A.Tonazzo 9 , L.Tortora 38 , G.Transtromer 24 , D.Treille 9 , G.Tristram 8 , M.Trochimczuk 51 , C.Troncon 27 



iii 



A.Tsirou 9 , M-L.Turluer 39 , I.A.Tyapkin 16 , S.Tzamarias 11 , B.Ueberschaer 52 , O.Ullaland 9 , V.Uvarov 42 , G.Valenti 5 
E.Vallazza 46 , G.W.Van Apeldoorn 30 , P. Van Dam 30 , W.K.Van Doninck 2 , J. Van Eldik 30 , A. Van Lysebetten 2 
I. Van Vulpen 30 , N.Vassilopoulos 34 , G.Vegni 27 , L.Ventura 35 , W.Venus 36 ' 9 , F.Verbeure 2 , M.Verlato 35 , L.S.Vertogradov 16 
V.Verzi 37 , D.Vilanova 39 , L.Vitale 46 , EVlasov 42 , A.S.Vodopyanov 16 , C.Vollmcr 17 , G.Voulgaris 3 , V.Vrba 12 , H.Wahlen 52 
C.Walck 44 , C.Weiser 17 , D.Wicke 52 , J.H.Wickcns 2 , G.R.Wilkinson 9 , M.Winter 10 , M.Witek 18 , G.Wolf 9 , J.Yi 1 
O.Yushchcnko 42 , A.Zaitsev 42 , A.Zalewska 18 , P.Zalewski 51 , D.Zavrtanik 43 , E.Zevgolatakos 11 , N.I.Zimin 16 ' 24 
G.C.Zucchelli 44 , G.Zumerle 35 



1 Department of Physics and Astronomy, Iowa State University, Ames IA 50011-3160, USA 
2 Physics Department, Univ. Instclling Antwerpen, Universiteitsplein 1, BE-2610 Wilrijk, Belgium 
and IIHE, ULB-VUB, Pleinlaan 2, BE-1050 Brussels, Belgium 

and Faculte des Sciences, Univ. dc l'Etat Mons, Av. Maistriau 19, BE- 7000 Mons, Belgium 
3 Physics Laboratory, University of Athens, Solonos Str. 104, GR-10680 Athens, Greece 
4 Department of Physics, University of Bergen, Allegaten 55, NO-5007 Bergen, Norway 
5 Dipartimento di Fisica, Universita di Bologna and INFN, Via Irnerio 46, IT-40126 Bologna, Italy 
6 Centro Brasileiro dc Pcsquisas Fi'sicas, rua Xavier Sigaud 150, BR-22290 Rio de Janeiro, Brazil 
and Depto. de Fisica, Pont. Univ. Catolica, CP. 38071 BR-22453 Rio de Janeiro, Brazil 
and Inst, de Fisica, Univ. Estadual do Rio de Janeiro, rua Sao Francisco Xavier 524, Rio de Janeiro, Brazil 
7 Comenius University, Faculty of Mathematics and Physics, Mlynska Dolina, SK-84215 Bratislava, Slovakia 
8 College de France, Lab. de Physique Corpusculaire, IN2P3-CNRS, FR-75231 Paris Cedex 05, France 
9 CERN, CH-1211 Geneva 23, Switzerland 

10 Institut de Recherches Subatomiques, IN2P3 - CNRS/ULP - BP20, FR-67037 Strasbourg Ccdcx, France 
"Institute of Nuclear Physics, N.C.S.R. Dcmokritos, P.O. Box 60228, GR-15310 Athens, Greece 

12 FZU, Inst, of Phys. of the C.A.S. High Energy Physics Division, Na Slovancc 2, CZ-180 40, Praha 8, Czech Republic 
13 Dipartimento di Fisica, Universita di Genova and INFN, Via Dodecaneso 33, IT-16146 Genova, Italy 
14 Institut des Sciences Nucleaires, IN2P3-CNRS, Universite de Grenoble 1, FR-38026 Grenoble Cedex, France 
15 Helsinki Institute of Physics, HIP, P.O. Box 9, FI-00014 Helsinki, Finland 

16 Joint Institute for Nuclear Research, Dubna, Head Post Office, P.O. Box 79, RU-101 000 Moscow, Russian Federation 

17 Institut fur Expcrimcntellc Kcrnphysik, Univcrsitat Karlsruhe, Postfach 6980, DE-76128 Karlsruhe, Germany 

18 Institute of Nuclear Physics and University of Mining and Mctalurgy, Ul. Kawiory 26a, PL-30055 Krakow, Poland 

19 Universite de Paris-Sud, Lab. de l'Accelerateur Lineaire, IN2P3-CNRS, Bat. 200, FR-91405 Orsay Cedex, France 

20 School of Physics and Chemistry, University of Lancaster, Lancaster LAI 4YB, UK 

21 LIP, 1ST, FCUL - Av. Elias Garcia, 14-1°, PT-1000 Lisboa Codex, Portugal 

22 Department of Physics, University of Liverpool, P.O. Box 147, Liverpool L69 3BX, UK 

23 LPNHE, IN2P3-CNRS, Univ. Paris VI et VII, Tour 33 (RdC), 4 place Jussieu, FR-75252 Paris Cedex 05, France 

24 Department of Physics, University of Lund, Solvcgatan 14, SE-223 63 Lund, Sweden 

25 Universite Claude Bernard dc Lyon, IPNL, IN2P3-CNRS, FR-69622 Villeurbanne Cedex, France 

26 Univ. d'Aix - Marseille II - CPP, IN2P3-CNRS, FR-13288 Marseille Cedex 09, France 

27 Dipartimento di Fisica, Universita di Milano and INFN, Via Celoria 16, IT-20133 Milan, Italy 

28 Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

29 NC, Nuclear Centre of MFF, Charles University, Areal MFF, V Holcsovickach 2, CZ-180 00, Praha 8, Czech Republic 
30 NIKHEF, Postbus 41882, NL-1009 DB Amsterdam, The Netherlands 

31 National Technical University, Physics Department, Zografou Campus, GR-15773 Athens, Greece 
32 Physics Department, University of Oslo, Blindern, NO-1000 Oslo 3, Norway 
33 Dpto. Fisica, Univ. Oviedo, Avda. Calvo Sotclo s/n, ES-33007 Oviedo, Spain 
34 Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK 
35 Dipartimento di Fisica, Universita di Padova and INFN, Via Marzolo 8, IT-35131 Padua, Italy 
36 Rutherford Appleton Laboratory, Chilton, Didcot OX11 OQX, UK 

37 Dipartimento di Fisica, Universita di Roma II and INFN, Tor Vergata, IT-00173 Rome, Italy 

38 Dipartimento di Fisica, Universita di Roma III and INFN, Via della Vasca Navale 84, IT-00146 Rome, Italy 

39 DAPNIA/Service de Physique des Particules, CEA-Saclay, FR-91191 Gif-sur-Yvette Cedex, France 

4(1 Instituto de Fisica de Cantabria (CSIC-UC), Avda. los Castros s/n, ES-39006 Santander, Spain 

41 Dipartimcnto di Fisica, Universita degli Studi di Roma La Sapienza, Piazzalc Aldo Moro 2, IT-00185 Rome, Italy 

42 Inst. for High Energy Physics, Serpukov P.O. Box 35, Protvino, (Moscow Region), Russian Federation 

43 J. Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia and Laboratory for Astroparticle Physics, 

Nova Gorica Polytechnic, Kostanjcviska 16a, SI-5000 Nova Gorica, Slovenia, 

and Department of Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 
44 Fysikum, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden 

45 Dipartimcnto di Fisica Spcrimcntale, Universita di Torino and INFN, Via P. Giuria 1, IT-10125 Turin, Italy 
46 Dipartimento di Fisica, Universita di Trieste and INFN, Via A. Valerio 2, IT-34127 Trieste, Italy 

and Istituto di Fisica, Universita di Udine, IT-33100 Udine, Italy 
47 Univ. Federal do Rio dc Janeiro, CP. 68528 Cidade Univ., Ilha do Fundao BR-21945-970 Rio dc Janeiro, Brazil 
48 Department of Radiation Sciences, University of Uppsala, P.O. Box 535, SE-751 21 Uppsala, Sweden 
49 IFIC, Valencia-CSIC, and D.F.A.M.N., U. de Valencia, Avda. Dr. Molincr 50, ES-46100 Burjassot (Valencia), Spain 
50 Institut fur Hochcncrgicphysik, Osterr. Akad. d. Wisscnsch., Nikolsdorfcrgasse 18, AT-1050 Vienna, Austria 
51 Inst. Nuclear Studies and University of Warsaw, Ul. Hoza 69, PL-00681 Warsaw, Poland 
52 Fachbereich Physik, University of Wuppertal, Postfach 100 127, DE-42097 Wuppertal, Germany 

53 On leave of absence from IHEP Serpukhov 

54 Now at University of Florida 



1 



1 Introduction 

The cross-section for the process e + e~ — > Z — > // for the fermion / as a function of 
its polar angle 9 with respect to the direction of the e~ can be expressed as: 

d°~ 8 a f n In 

oc 1 + - A J FB cos 9 + cos 9 . 



d cos 9 3 

The term proportional to cos 9 generates a forward-backward asymmetry A^ B which re- 
sults from the interference of the vector (v) and axial vector (a) couplings of the initial 
and final state fermions of the Z boson. The improved Born level asymmetry using the 
effective couplings for pure Z exchange is given by: 

f 3 2v e a e 2vf(if 

FB = l ^Ts 2 ^Ts 2 

The measurement of the forward-backward asymmetry for the different fermions in the 
final state can thus be used to measure v/a, and hence to determine the electroweak 
mixing angle sin 2 #e//- 

In this analysis, the forward-backward asymmetries for the processes e + e~ — > cc and 
e + e~ — > bb at the Z resonance are measured using reconstructed D mesons in the modes 1 : 

D*+ -> D°n + 

-> {K-it + )it + 

— > (K~ n + n~ n + )n + 

— > (K^7i + (7i ))n + with and without n° reconstruction 
-> (X-// + (^))tt + 
-> (K-e + (is e ))n+ 

D° — > i^ _ 7T + (7r°) without 7T° reconstruction 

The D meson contains a charm quark and therefore provides a clean signature of a cc 
event or a decay of a heavy 6-hadron in a 66 event. In both cases the charge state 2 of the 
D is directly correlated to the charge of the primary quark. 

Particle identification in the DELPHI detector is provided by ring imaging Cherenkov 
counters (RICH), the specific energy loss dE/dx in the Time Projection Chamber (TPC), 
and electron and muon identification (see section 2); together with 7r° reconstruction, it 
is used to identify the D decay products and to reduce the combinatorial background. 

It is necessary to distinguish between the contributions of D mesons from c and b 
quark events in order to determine the individual forward-backward asymmetries. In 
this analysis a simultaneous fit is performed in terms of the scaled energy Xe = 2Ed/ ' yfs 
(where Ed is the energy of the D) and the 6-tagging probability variable (see section 5). 
The forward-backward asymmetry is then extracted from the distribution of the cosine 
of the polar angle of the thrust axis signed by the charge state of the D meson. 

In this paper an update of previous DELPHI results [1] is presented. All LEP 1 data 
collected by the DELPHI detector in the years 1992 to 1995 are used. The high quality 
of the data through the final LEP 1 reprocessing in combination with improved recon- 
struction and 6-tagging techniques [2] result in a significant gain in statistical precision. 

1 Throughout the paper charge-conjugated states are included implicitly. 

2 For the D° the charge state is defined as minus the charge of its decay kaon. 
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2 The DELPHI Detector 

The DELPHI detector consists of several independent devices for tracking, calorimetry 
and lepton and hadron identification. The components relevant for this analysis will be 
briefly described in the following. A detailed description of the whole apparatus and its 
performance can be found in [3]. 

In the barrel region the innermost component is the vertex detector (VD) near to the 
LEP beam pipe. The VD consists of three concentric layers (closer, inner and outer) 
of silicon microstrip detectors. Since 1994 the VD provides R<p and z information 3 in 
the closer and outer layer and has an extended polar angle coverage of 25° < 9 < 155° 
(closer layer). For polar angles of 44° < 9 < 136°, a particle crosses all three layers of 
the VD. With an intrinsic R(f> resolution of 7.6 fim [3], the VD is the main component for 
reconstructing secondary vertices of heavy hadron decays. 

The inner detector (ID) is outside the VD and consisted of a jet-chamber to perform a 
precise R<p measurement, and five cylindrical MWPC layers. In 1995 the MWPC layers 
were replaced by five layers of straw tube detectors. 

The ID is followed by the TPC, the main tracking device in DELPHI. It covers polar 
angles between 21° < 9 < 159° with a single point resolution for charged particles of 
approximately 250 /im in Rep and 880 /im in z [3] . The analysis of the pulse heights of 
the signals of up to 192 sense wires allows the determination of the specific energy loss, 
dE/dx, of charged particles which can be used for particle identification (see section 4). 

The barrel ring imaging Cherenkov counter is behind the TPC. Its gas and liquid 
radiators allow particle identification for pions, kaons and protons over almost the whole 
momentum range (see section 4). 

The outer detector (OD) is mounted behind the RICH to give additional tracking 
information. It improves significantly the momentum resolution due to its large distance 
from the interaction point. Five layers of drift cells cover polar angles between 42° < 9 < 
138° and provide R(j> and z information. 

The barrel electromagnetic calorimeter (HPC) is between the OD and the supercon- 
ducting coil and covers polar angles between 42° < 9 < 138°. It is a gas-sampling device 
which provides complete three-dimensional charge information in the same way as a time 
projection chamber. The excellent granularity allows good separation between close par- 
ticles in three dimensions. This permits good electron identification even inside jets and 
direct identification of tt° — > 77 decays. 

In the forward region, tracking is performed by two planar drift chambers (FCA and 
FCB) with a polar angle covering of 11° < 9 < 33° (FCA) and 11° < 9 < 36.5° (FCB). 
Their resolutions transverse to the beam axis are 270 fim (FCA) and 150 /im (FCB) 
respectively. 

For muon identification the DELPHI detector is surrounded by layers of drift chambers. 
They cover 52° < 9 < 128° with a resolution of 4 mm in R<p and 2.5 cm in z in the barrel 
region and 9° < 9 < 43° with a resolution of 1 mm in the forward region. 

3 Hadronic event selection 

Charged particles are selected as follows. The momentum is required to be between 
0.4 GeV/c and 50 GeV/c, the relative error on the momentum measurement less than 
100%, the polar angle relative to the beam axis between 20° and 160°, the length of 

3 In the DELPHI coordinate system, z is along the electron beam direction, (j> and R are the azimuthal angle and radius 
in the xy plane, and 8 is the polar angle with respect to the z axis. 
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tracks in the TPC larger than 30 cm, the projection of the impact parameter relative to 
the interaction point less than 4 cm in the plane transverse to the beam direction, and 
the distance along the beam direction to the interaction point less than 10 cm. 

Hadronic events are selected by requiring five or more charged particles and a total 
energy in charged particles larger than 12 % of the collision energy (assuming all charged 
particles to be pions). A total of 3.5 million hadronic events is obtained from the 1992- 
1995 data, at centre-of-mass energies within ± 2 GeV of the Z resonance mass. According 
to the simulation, the selection efficiency for hadronic Z decays is 95.7%. Table 1 shows 
the number of selected hadronic events. Remaining backgrounds from r pairs and Bhabha 
events are found to be negligible for this analysis. 

A set of about 8.5 million simulated hadronic events for the years 1992 to 1995 is used. 
They are generated using JETSET 7.4 Parton Shower model [4] in combination with the 
full simulation of the DELPHI detector. The parameters of the generator are tuned to 
the DELPHI data [5]. 

For each event, the primary interaction vertex is determined from the measured tracks, 
with a constraint from the measured mean beam spot position. The removal of the track 
with the largest \ 2 (followed by a refit of the vertex) is repeated until either the \ 2 °f 
each contributing track is less than 3 or less than three charged particle tracks are left. 
All track parameters are recalculated after a helix extrapolation to this vertex position. 
The resolution of tracks measured only by the forward tracking chambers is improved by 
a track refit using the primary vertex. Forward tracks having a x 2 in the refit larger than 
100 are removed from the analysis. 



Year 


91.235 GeV 


Data events 
89.434 GeV 


92.990 GeV 


Simulation 

events 


1992 


703859 






2003142 


1993 


475151 


97623 


134240 


1893139 


1994 


1386191 






3551362 


1995 


458700 


84763 


131637 


1126557 


92-95 


3023901 


182386 


265877 


8574200 



Table 1: The number of selected hadronic events for data and simulation. 



4 Reconstruction of charmed mesons 

Reconstructed D mesons are used as a signature for cc and bb events and are identified 
through their decay products. The D mesons are reconstructed in nine different decay 
modes (see table 2). In the following a brief description of the selection criteria for the 
D candidates is given. 

For all decay modes the selection of candidates is performed in a similar way. A 
number of charged particles (corresponding to the multiplicity of the specific D°^ + decay 
mode) with momentum p larger than 1 GeV/c are combined, requiring the total charge 
to be zero in case of the D° and one in case of the D + decay. The invariant mass m D 
of the D°/ + candidate is calculated, assuming one of the particles to be a kaon and the 
others pions. In addition the kaon momentum has to exceed 2 GeV/c for the D + decay, 
the leptonic modes and the decays with 7r° reconstruction, and 1 GeV/c for all the other 
decay channels. A D* + candidate is obtained by associating a low momentum pion down 



4 



to 0.4 GeV/c to the reconstructed D° meson. The charge of the pion is required to be 
opposite to that of the kaon from the D° decay. 

For the semileptonic decay modes D° — > K~e + v and D° — > K~/j, + v, the lepton is 
required to be identified, using standard DELPHI identification criteria [3,6]. 

For the reconstruction of 7r° — > 77 decays, three different classes of candidates mea- 
sured in the HPC are used [3,7]. The measurement of two separated photon showers in 
the HPC are used to construct the 7T°. Above n° energies of 6-8 GeV, the angle between 
the 77 pair is too small to separate the photon showers. The pion is therefore derived 
from the analysis of the shower shape of these merged photons. Information from photons 
converted in front of the TPC is also used to reconstruct photons and thereby the neutral 
pion. 

The particle identification provided by the RICH and the specific energy loss dE/dx 
measurement in the TPC are used to reduce the combinatorial background. Due to 
the large number of pions in the hadronic final state, combinations in which a pion is 
assigned as a kaon candidate are the main contribution to the background. To optimize 
the efficiency of the D signal, a pion veto, rather than direct kaon identification, has 
been introduced. Tagging is performed using DELPHI standard tagging routines for 
the RICH [8] and the dE/dx [3] identification. For the RICH, the measured Cherenkov 
angle information is translated into ir, K and p tagging information tagpj;^ p , taking into 
account the quality of the measurement. The way the information from the two radiators 
is combined depends on the momentum of the candidate, in order to guarantee the best 
separation over almost the whole momentum range. Kaon candidates are tagged if they 
have either no pion tag or a very loose pion tag tag^ ich . 

The dE/dx information is used only if no RICH information is available. For each 
track a probability P can be expressed in terms of the expected ionisation for a given 
particle hypothesis: 

Pk„ = exp^-\ ^E/dx-dE/dx^ 2 ) ^ 

This can be translated into a normalized kaon probability: 

p 

tag TPC = * (2) 

Depending on the decay channel, a cut tag TPC > 0.2 — 0.3 to the kaon candidate is 
applied. 

For the D° and the D + decays the kaon has to be tagged by the RICH or the TPC 
as explained above. For the D* + decay channels a tag of the kaon candidate is not 
necessarily required because of the cut on the mass difference Am = m D *+ — m D o which 
selects rather pure sample of D mesons. 

For all decay modes a secondary vertex fit for the D°^ + is performed and the D°^ + 
flight distance and improved track parameters are obtained. All tracks associated to a 
D are required to have at least one hit in the vertex detector. A further reduction of 
background for the D + is achieved by rejecting track combinations with a probability for 
X 2 of the vertex fit less than 0.001. The slow pion from the D* + decay is constrained to 
the D° vertex, which is a good approximation for the D* + decay vertex because of the 
small transverse momentum of the slow pion with respect to the flight direction of the 
D°. 

Cuts on the helicity angle distribution are used to achieve a further significant reduc- 
tion of the combinatorial background. The helicity angle 9h is defined as the angle of 
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decay mode 


AL for D 


cm 




E 


a 


b 


c 


d 


D*+ -> (K-7i + )n + 


-0.05 


to 2 


.0 


0.15 


-2.5 


0.04 


3.0 


0.1 


D*+ -> (K-7r+7r-7r+)7r+ 


0.0 


to 2 


.0 


0.3 


-1.0 


0.05 


2.0 


0.1 




0.0 


to 2 


.0 


0.3 


-1.5 


0.09 


2.0 


0.1 


-> (K-fi + iy)TT + 


0.0 


to 2 


.0 


0.2 


-1.2 


0.06 


3.0 


0.1 


D*+ -> (^-e + z/)7r+ 


0.0 


to 2 


.0 


0.2 


-1.4 


0.05 


3.0 


0.1 


D*+ -> (^-7T+(7r ))7T + 


0.0 


to 2 


.0 


0.2 


-2.5 


0.03 


3.0 


0.0 


D+ ^^-7T + 7T + 


0.125 


to 2 


.0 


0.35 






3.0 


0.1 


L>° -> + 


0.05 


to 2 


.0 


0.3 


-0.5 


0.125 


2.0 


0.2 


D° ^^-7T+(7T ) 


0.05 


to 2 


.0 


0.3 


-0.6 


0.15 


2.0 


0.2 



Table 2: The selection cuts for X E versus AL (a,b), X E versus cos9 H (c,d) and the 
absolute cut on X E - 

the sphericity axis in the L> 0//+ rest frame with respect to the L> 0//+ flight direction. The 
orientation of the sphericity axis corresponds to the kaon candidate in the decay. L>°/ + 
decays are isotropic in cos9 H , whereas the background is peaked at cos6 H = ±1- Be- 
cause of the shape of the energy spectrum of charged particles in hadronic Z events, the 
combinatorial background is concentrated at small scaled energies Xe(D). Therefore X E 
cuts which depend on the helicity angle are used, in order to remove higher background 
contributions at small D meson energies: 

X E > 0.5 • e 6 *' 00 ***'- 1 ) +d. (3) 

In addition the scaled energy X E of the D combination is required to exceed the limits 
Xf* given in table 2, where the parameters c and d are also shown. 

The distance between the primary and the L>°/+ vertex is calculated in the plane 
transverse to the beam axis and projected onto the D°/ + direction of flight to obtain the 
decay length AL. A vertex combination is accepted if AL is within the range specified in 
table 2 for the different decay modes. In addition an X E dependent cut on AL is applied: 

AL(X E )>a-(X E -X c ^) 2 + b. (4) 

This rejects combinatorial background, which is concentrated at low values of both of 
these variables. The parameters a and b are given in table 2. 

For the D + — > K~n + n + mode, a cut 4 of Am > 200 Mev is used to veto D* + decays. 

For the D* + — > (7^~7r + 7r°)7r + decay mode with a reconstructed 7r°, an additional cut 
on the D° Dalitz plot is applied to select the dominant decay via D° — > K~p + : 

0.5 GeV/c 2 < m K - n+ < 1.1 GeV/c 2 and 1.4 GeV/c 2 < m K - n o < 1.8 GeV/c 2 
1.4 GeV/c 2 < m K - n+ < 1.8 GeV/c 2 and 0.5 GeV/c 2 < m K - w o < 1.1 GeV/c 2 

The mass bands to select the different D°/ + decay modes and the cuts on the mass 
difference are listed in table 3. The mass and mass difference bands in the signal regions 
for the different decay channels used in the final analysis are given in table 4; the sidebands 
are defined as the mass or mass difference regions given in table 3 not selected by the cuts 

4 The mass difference for the D + is given as Am = m D+ — m K - lr +, in analogy with that for the D* decay: Am = 
ni£>» — m D o. Both possible m j4 -- Jr + combinations are tested. 
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mode 



D o/+ 

mass interval 
[GeV/c 2 ] 



max. Am 
[GeV/c 2 ] 



D*+ -> (^-7r+)7r+ 

+ -> (^-7T + 7T-7r + )^ 
+ -> (ir-7T + 77)7T + 

-> (K-^ + v)-n + 
D*+ -> (^-e+z/)7r+ 
D*+ -> (ir-vr + (7r ))7r + 

D+ ^^'7T + 7r + 
D° ^^-7T+(7T ) 



1 


790 


to 


1 


94 


1 


845 


to 


1 


90 


1 


740 


to 


1 


98 





750 


to 


1 


75 





750 


to 


1 


75 


1 


350 


to 


1 


75 


1 


700 


to 


2 


05 


1 


750 


to 


2 


20 


1 


500 


to 


1 


70 



0.160 
0.160 
0.165 
0.250 
0.250 
0.175 



Table 3: Mass and mass difference cuts for the selection of the D meson signal plus 
sideband regions. 



in table 4. The mass differences for the D* + and also the D° and D + mass distributions 
are shown in figures 1 and 2. The wrong sign combination for D° — > K~tt + , where the 
kaon and pion masses are interchanged, is also shown in figure 2. The histograms show 
the simulated distributions normalized to the data samples. The contributions of signal 
and background are adjusted to compensate for different D rates in data and simulation. 



5 Measurement of A C FB and A h FB 

For a measurement of A C FB and A b FB from the polar angle cos 6 t hrust of the thrust axis 
in the D meson events, it is necessary to separate D from cc and bb events and the com- 
binatorial background. Since the c and b asymmetries are expected to be of comparable 
size and to have the same relative sign, the statistical precision of the measurement is 
limited by the negative correlation between both asymmetries. In this analysis, good 
separation with a small correlation is obtained by using the scaled energy distribution 
X E of the D candidates and the event 6-tagging variable V ev [2]. 

The hadronization of primary c quarks leads to high energy D mesons, whereas b 
quarks fragment into 6-hadrons which then decay into D mesons with a softer energy 
spectrum. The combinatorial background is concentrated at low Xe- Furthermore bb 
events can be identified by 6-tagging, which utilises special features of B hadrons, as 
compared with other hadrons. The combined 6-tagging used in this analysis takes into 
account the long lifetime and the large mass of B hadrons, their higher decay multiplicity 
and their large Xe(B). This leads to a high tagging efficiency in combination with good 
separation power of the tag. 

The shape of the combinatorial background is tested using the sidebands in the mass 
(or the mass difference) distribution. Due to the different relative acceptance of D mesons 
and background at small and large polar angles, the fit method has to take into account 
the | cos Othrust | dependence of the different classes. 

The charge state of a signal D is directly correlated to the charge of the primary quark, 
whereas the charge correlation of the combinatorial background is expected to be very 
small. 
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For the asymmetry measurement, partially reconstructed D* + mesons {tt s i + X) and 
reflections from other decay modes (see figures 1 and 2) have to be considered as signal to 
avoid charge correlations in the background. The contributions from reflections, where 
some particles from the D decay are assigned a wrong mass or are missing, and true 
D decays are treated as one class, because of the similar shape of the signals and the 
charge correlation with the primary quark. This leads to a significant increase of the 
sample for the K~-K + {it ) decay mode. The rate of partially reconstructed D* + mesons, 
where a tt + from a D* + decay is combined with a fake D°, depends on the branching 
ratio D* + — > D°ir + , the D* + production rate and the efficiency in the relevant mass 
difference interval. The contribution of partially reconstructed D* + decays to the signal 
is taken from the simulation and contributes to the systematic uncertainty. In the case 
of D° — > K~ir + /K~ir + (-K°) decay modes without the D* + constraint, candidates with 
wrong mass assignments flip the sign of the estimated primary quark direction. This is 
taken into account in the fit; the systematic error allows for uncertainties. 

To avoid double counting of events, only one D candidate in the signal region per 
event is retained. For a given event the D candidate with the largest kaon momentum is 
used. If two candidates for a given decay mode use the same kaon track, the one with the 
largest Xe{D) is used. Events entering the signal region for the K~n + decay mode are 
removed from the fT~7r + 7r~7r + distribution and events from both decay modes are then 
removed from the K~fi + v or K~e + v distribution and so forth. The order of the rejection 
is as listed in table 3. Good agreement between data and simulation was found for the 
rejection. 

The numbers of reconstructed D decays given in table 4 are obtained from fits to the 
mass spectra. A total sample of 61829 ± 521 reconstructed D decays is used for the 
asymmetry measurement. The D mass bands to select D meson candidates are listed in 
table 4. 



decay mode 



signal events 



signal region 
Am [GeV/c 2 ] 



signal region 
m D [GeV/c 2 ] 



Rs/B 


0.95 


± 


0.02 


0.86 


± 


0.02 


1.19 


± 


0.03 


0.64 


± 


0.02 


0.98 


± 


0.04 


1.16 


± 


0.02 


0.83 


± 


0.02 


1.00 


± 


0.02 


1.21 


± 


0.04 



D* + 

D*+ 
D *+ 

D *+ 

D *+ 

D *+ 

D+ 
D° 
D° 



(K-7i + )n + 

(K-/2 + u)7l + 

(K-e + is)ir + 
(K-n + (7i°))n + 

K-TI + TI + 
K~7T + 



6030 
5123 
5787 
3042 
1810 
15111 
5667 
9311 
9948 



± 103 
± 103 
± 125 
± 91 
± 65 
± 232 
± 161 
± 232 
± 298 



0.143-0.148 
0.143-0.148 
0.141-0.151 

< 0.180 

< 0.180 

< 0.152 



1.83-1.91 
1.80-1.93 
1.50-1.70 



Table 4: D meson samples used for the measurement, cuts to select signal regions, and 
the relative normalizations Rs/b of signal to background for data and simulation. 



5.1 The minimum % 2 fit 

The determination of the asymmetries at s/s = 91.235 GeV is achieved by a minimum 
X 2 fit to the D samples using the scaled energy X E) the transformed 6-tagging variable 
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Decay mode 


number of bins per 


average number 


v/s = 91.235 GeV 


X E 


tr(P e „) 


COS 9 thrust 


of events 


D*+ -> (K-n + )n + 


4 


5 


4 


77.8 


D*+ -> (K-7r + 7r-7r + )7r+ 


5 


5 


5 


65.8 




4 


4 


4 


88.0 


£>*+ -> (K-fi + is)ir + 


4 


6 


4 


49.6 


D*+ -> (K-e+z/)7r+ 


4 


4 


3 


56.4 


D*+ -> (K-7T + (7T ))7r + 


6 


7 


5 


77.7 


L> + ^^-7T + 7T + 


5 


5 


5 


87.9 




5 


6 


5 


80.0 


-D° ^K-7T+(7T ) 


7 


7 


7 


59.9 



Table 5: Number of bins in each dimension used for the individual decay modes and the 
average number of data events per bin at y/s = 91.235 GeV. 

tr(V e v) for the event and the polar angle Q-cos Q thrust signed with the charge state Q of the 
D. Examples of these distributions for the D* + — > (K~7t + )tt + channel are shown in figure 
3. The measured distributions are compared to the predictions of the simulation, split 
into charm, bottom and background events. The simulated prediction is normalized to 
the data to reproduce the signal to background ratio. Therefore a factor Rs/b (see table 
4) is introduced for each decay mode, which compensates for different D rates in data 
and simulation. After this correction, good agreement is found in all distributions. The 
shape of the background distribution, as obtained from the sidebands, is well reproduced 
by the simulation. 

A transformation of the event variable P ev is used for the 6-tagging distribution: 



The bins in the three dimensional Xe, tr(V ev ) and cos Qthrust space have been chosen such 
that each bin contains about 70 events (table 5). In each bin % the differential asymmetry: 

is calculated from the numbers of events and N~ with Q ■ cos Qthrust greater or less 
than zero, respectively. The observed asymmetry receives contributions from c, b and 
combinatorial background. The fractions fji of D signal and reflection from c and b 
events as well as the fractions of partially reconstructed D mesons and combinatorial 
background are taken from the simulation. Furthermore the combinatorial background 
is divided into c, b and uds contributions to account for the small charge correlation 
to the primary quark in the background, especially for the semileptonic decay mode of 
the D°. This leads to three different contributions (i.e. signal plus reflections, partially 
reconstructed D, and combinatorial background) from each of c and b, and one for the 
background from uds. 

The x 2 to be minimized is given by: 
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Nbins [ 7 

X 2 = £ A*'* - £ fa Cji A j FB (cos 9,) } ja] 



(7) 



i=i 



j'=i 



where <7j accounts for the statistical error of both data and simulation; A FB (cos9i) is the 
differential asymmetry: 



A^ B (cos Oi, 



A\ 



COS 6,; 



FB 



1 + COS 2 Oi ' 



(8) 



of 6, c or wds events; and Cji is the charge correlation of the class j calculated in each 
bin using the simulation. For b events the mixing effect leads to values of the charge 
correlation Cji which are smaller than 1, and thus to a smaller observed b asymmetry. 
The simulation is used to estimate the mixing effect as a function of the 6-tagging tr(P ev ), 
because the 6-tagging depends on the individual B lifetime (see section 6 for details). 
The combinatorial background from c is expected to have only a small charge correlation 
Cji to the primary quark at large energies X F - The asymmetry of the combinatorial 
background and of D mesons from gluon splitting in uds events is expected to be very 
small. The predictions from the simulation for this class in each bin are subtracted in the 
fit. The two fit parameters A C FB and A b FB are used for all three classes from c and b, while 
for the single class from uds the prediction of the simulation is used. The agreement of 
data and simulation is tested in the sidebands of the different samples where no significant 
deviations are found. 



Decay 



mode 
$9,434 GeV 



number of bins per 



X, 



tr(P e 



cosU thrust \ 



number of 
candidates 



D* + -> {K-tt + )tt + 
D*+ -> (K-it + it-it + )it- 
D*+ -> (fr-7r+77)7r+ 
D*+ -> (K-fi + is)ir + 
D*+ -> (K-e+is)ir+ 
D*+ -> (^-7T+(7r°))7r+ 

D+ -> K-1T + 1T + 



3 
4 
3 
3 
3 
5 
4 
4 
4 



3 
4 
3 
3 
3 
5 
4 
5 
5 



4 
3 
4 
4 
2 
4 
4 
3 
5 



378 
539 
356 
260 
188 

1096 
814 
721 

1287 



Table 6: Number of bins in each dimension used for the individual decay modes and the 
number of data events per bin at y/s = 89.434 GeV. 



5.2 The maximum likelihood fit 

The data taken by the DELPHI detector in the years 1993 and 1995 at energies 
near the Z resonance allow the investigation of the energy dependence of the forward- 
backward asymmetry of c and b quarks. Due to the reduced statistics of 182386 events at 
x/s = 89.434 GeV and 265877 events at y/s = 92.990 GeV, a binned maximum likelihood 
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Decay mode 


number of bins per 


number of 


y/s = 92.990 GeV 


X E 


tr(V ev ) 


COS 9 thrust 


candidates 


D*+ -> (K-7i+)n + 


3 


4 


4 


576 


D*+ -> (K-7r + 7r-7r + )7r+ 


4 


4 


4 


734 


-> (fs:-7r + 77)7r + 


3 


4 


3 


473 


-> (K-h+is)it + 


3 


3 


4 


412 


D*+ -> (K"e + z/)7r+ 


3 


3 


3 


276 


D*+ -> (K-7T + (7r ))7T + 


5 


5 


5 


1475 




5 


4 


4 


1058 




4 


5 


4 


1076 


D° ^K-7T+(7T ) 


5 


5 


6 


1947 



Table 7: Number of bins in each dimension used for the individual decay modes and the 
number of data events per bin at y/s = 92.990 GeV. 



fit is used instead of the \ 2 fit °f section 5.1. The likelihood function is given by: 



E 

i=i 



In 



A 



, NT 

+ » 



+ ln- K 



(9) 



where the Nj,i ns cells are given by the binned information in Xe, tr(V ev ) and cosOthmst- 
The bins for the three dimensions are given in tables 6 and 7. The A^ describe the 
expectation in each cell of a Poisson distribution, and are given by: 



A 



± 



fjtot 7 

-^-^/^■(l±^ B (cos^))- 



(10) 



The total number of candidates N, 



tot 



+ N; is taken from real data and the coeffi- 



cients fij, Cij and Ap B (cos9i) are defined in the same way as in section 5.1. 



6 Effective mixing in b — > D decays 

The observed forward-backward asymmetry in b events is proportional to the charge 
correlation of the reconstructed D meson to the primary quark. The correlation C& is 
reduced because of two factors. The first is B° — B° mixing, and the second is double D 
production in B decays. The observed b asymmetry in the different D samples needs to 
be corrected for both of the above sources of D mesons of the wrong charge state. 

The effect of B° — B° mixing is to reduce the correlation C b by a factor 1 — 2%. The 
mixing probability x is determined by the mass difference Am between the two mass 
eigenstates and by the B lifetime tb- The product of these leads to Xd — 0.172 ± 0.010 
[9]. For the B° s only a lower limit of Am s > 14 ps -1 [9] is known. This is compatible with 
full mixing Xs ~ 0.5. 

The production of D mesons from the "upper vertex" (i.e. via the W decay, rather 
than from the bcW coupling) also reduces the charge correlation to the primary b quark. 
A sizeable rate J'b^w^d of wrong sign D mesons from B decays reduces the measured 
b asymmetry by a factor 1 — 2 • Jb^w^d- Recent measurements of CLEO and ALEPH 
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indicate a significant rate of double-charmed B decays involving no D s production (D s 
production would not result in D or D* of the wrong charge state). CLEO [10] finds for 
a mixture of and B + a ratio of: 

^—^1 = 0.100 ± 0.026 (stat) ± 0.016 (syst). (11) 
T[B — ► DX) 

This number is in good agreement with the ALEPH result on double D production in B 
decays [11]: 

BR(b -> D°D°, D°D~, D + D° + X) = (7.8 1?;° t\l ^ j) %, (12) 

where the dominant contribution is B — > D^D^K. The first error is statistical, the 
second one contains the systematic errors and the third corresponds to the uncertainties 
on the D branching fractions. Both measurements include Cabbibo suppressed W — > cd 
decays which are expected to contribute with a rate of about 1% per W decay to the 
"upper vertex" charm rate. 

Mixing is relevant for D mesons from B® and B® decays, but not from B + or A& 
decays, whereas the effect from the "upper vertex" charm contributes in all B decays. 
The fractions of D ± , D° and D*^ from different B states need to be determined from 
the branching rates B — > D and B — > D to be able to consider the combination of both 
effects correctly. Very little is known at present about the individual exclusive branching 
ratios, but several inclusive measurements from the T(4S) and LEP experiments can be 
used to deduce the rates. 

CLEO and ARGUS_have measured the rates [9] of D°, D ± , Df and Af as well as 
the rates of D #± and 5*° in T(AS) decays, i.e. in decays of B% and B + at about 50% 
admixture. From these measurements the overall rates of B® and B + decays in the D 
samples are deduced, taking into account the production fractions of B + , 5°, B° s and A& 
in Z — > bb events [9]. The relative rate of D*^ from B~ and from B° is not measured. 
Therefore the JETSET prediction of 

(B + -> D*± X)/(B + >° -> D*± X) = 0.30 (13) 

with a relative error of 50 % is used. This number is compatible with the assumption 
that most of the B + — > D* ± X decays are produced via D** (i.e. D%, D{, D 1 , D ) decays 
and higher D resonances. 

The D* rates measured at the T(AS) also fix the effective rate of vector and pseu- 
doscalar mesons V/(V + P). The decay of the D*° into D + n~ is forbidden by phase 
space, while the branching ratio D* + — > D°tt + is measured to be 0.683 ± 0.013 [9]. This 
difference in charged and neutral D* decays significantly affects the fractions of B + and 
B® decays seen in the D ± and D° samples and thus the effective mixing. 

A small correction to the B + and B® into D ± , 5° and D*^ 1 rates originates from D** 
production. These states subsequently decay into vector or pseudoscalar charm mesons. 
D\ mesons decay into D* and D states and the ratio between these two decays is given 
by the measurement BR(D*° -> D + tt-)/BR(D*° -> D*+ti-) = 2.3 ± 0.8 [12]. Angular 
momentum conservation allows D\ and D\ mesons to decay into D* states, and D 
to decay into D states. The decay rates of the different D** states into charged and 
neutral states (D*7r,D7r) are fixed by isospin invariance. The relative production rates 
of the four D** states are assumed to be proportional to the number of spin states. The 
total D** rate is obtained from the measured semileptonic B branching ratios BR(B — > 
D**l+u) /BR(B -> Xl+u) = 0.26 ± 0.07 [9]. 
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The decay of B° s and A& into D ± 1 D° and D*^ also contribute to the sample. They 
can be deduced from the total rate of D°, D + , and A+ in Z — > bb events measured at 
LEP [13]. Here the number of charm quarks produced per b decay is limited to the LEP 
measurement n c = 1.17±0.04 [9], with charmonia and S c production taken into account. 

From this information and assuming that the branching fraction of W — > cs — > 
is equal to that for W — > cs — > the effective i? mixings in the D + , D° and -D* + 

samples are: 

Xeff(D + ) = 0.222 ±0.033 

X eff(D°) = 0.176 ±0.030 (14) 
X eff(D* + ) = 0.222 ± 0.033 . 

The errors quoted represent the precision of the measurements used to determine the b 
decay properties. The effective mixing for the D* + sample is in good agreement with a 
direct measurement of OPAL using the jet charge technique in the hemisphere opposite 
to the reconstructed meson. They obtained Xeff(D* + ) = 0.191 ± 0.083 [14]. 

7 QCD corrections 

The analysis of the final state of hadronic Z decays gives only indirect information 
about the electroweak process Z — > qq. The evolution to the final parton level and 
the following process of hadronization smear the clear signature of the initial qq system. 
The hard gluon radiation in the parton shower changes significantly the direction of 
the primary quarks and thus also the angular distribution of the following partons and 
hadrons. The size of this QCD effect strongly depends on the individual techniques in the 
determination of the forward-backward asymmetry. The QCD correction can be written 
as [15]: 

A FB = (1 - C q ) (A FB ) = (l-s q C£ CD ) (A FB ) (15) 

where (A fb )q is the asymmetry without gluon radiation, which can be calculated from 
the measured asymmetry A q FB through the correction coefficient C q . This correction 
coefficient can be parameterized by a bias factor s q , which accounts for the individual 
sensitivity to the QCD correction Cq CD . The values of the QCD corrections are estimated 
to be [15]: 

C^ CD = (2.96 ± 0.40)% 

Cqcd = (3-57 ± 0.76) % . (16) 

The experimental bias is studied using a fit to the simulation after setting the generated 
asymmetry to 75%. The observed relative difference C q of (—1.66 ± 0.01) % for cc and 
(-2.22 ± 0.02) % for bb lead to the bias factors s c = -46 % and s b = -75 %. The values 
of C q are used to define the experimental bias on the QCD corrections to correct the fit 
results. The statistical error of the fit and the uncertainty of the QCD correction are 
used to determine the systematic error. 

8 Systematic uncertainties 

The systematic error sources are of two types. The uncertainty of the simulation 
modelling of heavy quark production affects the measurement, and the fit method itself 
is a potential source for a systematic error. 
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systematic 




91.235 GeV 


89.434 GeV 


92.990 GeV 


error source 


variation 


SA C FB 


SA b FB 




SA b FB 


M C FB 


SA b FB 






xlO 3 


xlO 3 


xlO 3 


xlO 3 


xlO 3 


xlO 3 


MC statistics 


see text 


±2.47 


T3.54 


±3.31 


T7.20 


±3.31 


T7.20 


(X E ) D * 


0.510 ± 0.009 


±0.16 


T0.40 


tO.h 


T0.36 


±0.59 


T0.60 


(X e )b 


0.702 ± 0.008 


T0.20 


±0.25 


±0.69 


TO. 51 


T0.25 


T0.43 


e B^D 


0.42 ± 0.07 


±0.26 


T0.40 


T0.78 


T0.18 


±0.42 


±0.98 


t{B+) 


1.65 ± 0.04 


T0.06 


±0.18 


T0.02 


±0.31 


=fo.io 


±0.49 


r(B°) 


1.56 ± 0.04 


=F0.01 


T0.69 


±0.09 


T0.20 


^0. 11 


T0.28 


r(B° s ) 


1.54 ± 0.07 


=fo.oi 


T0.03 


±0.02 


±0.04 


T0.04 


±0.08 


r(A 6 ) 


1.22 ± 0.05 


T0.03 


±0.11 


±0.08 


T0.02 


T0.08 


±0.11 


t(D+) 


1.057 ± 0.015 


±0.02 


T0.03 


±0.01 


T0.09 


±0.02 


T0.03 


r(D°) 


0.415 ± 0.004 


±0.04 


±0.02 


T0.25 


±0.10 


±0.03 


±0.13 


T(Df) 


0.467 ± 0.017 


T0.02 


±0.04 


=fo.io 


±0.11 


T0.08 


±0.11 


r(A c ) 


0.206 ± 0.012 


±0.03 


T0.04 


T0.03 


±0.04 


T0.03 


±0.04 


f{D+) 


0.221 ± 0.020 


T0.02 


±0.04 


T0.35 


±0.46 


T0.27 


±0.36 


f(Df) 


0.112 ± 0.027 


±0.19 


T0.19 


T0.42 


±0.42 


±0.13 


T0.19 


f {Cbaryon) 


0.084 ± 0.022 


±0.03 


T0.03 


T0.23 


T0.05 


±0.20 


±0.15 


W-g^cc 


2.38 ± 0.48 % 


±0.05 


±0.19 


T0.05 


±0.19 


±0.05 


±0.19 


(R C -P C ^ D ) 


see text 


±0.38 


T0.40 




±0.07 


±0.70 


±0.22 


eff. mixing 


see text 


T0.05 


±5.78 


T0.26 


±2.82 


TO. 18 


±5.71 


QCD bias 


see text 


T0.24 


T0.23 


T0.17 


T0.18 


T0.41 


T0.27 


fit method 


see text 


±1.73 


T2.81 


±1.73 


T2.81 


±1.73 


T2.81 


hag 


see text 


±0.74 


±1.54 


±0.74 


±1.54 


±0.74 


±1.54 


Rs/B 


± 5/10% 


T0.63 


=Fl-05 


±1.81 


Tl.28 


T2.51 


T2.54 


n slow , wrong sign 


± 15% 


±0.21 


±0.43 


=Fl.07 


±0.97 


±0.48 


±0.85 


Af^ibackgr.) 


± 30% 


T0.55 


±0.02 


Tl.46 


T0.50 


T0.69 


T0.54 


A%(backgr.) 


± 30% 


Tl.26 


T3.62 


±2.04 


T7.86 


T3.93 


T6.59 


total 




±3.51 


T8.47 


±5.28 


Tll.67 


±6.16 


T12.20 



Table 8: Contributions to the systematic errors on the measured asymmetries. 
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To describe adequately heavy quark production, it is necessary to correct for inad- 
equate simulation settings. This is achieved using JETSET; the relevant distributions 
are compared (at the full simulation level, but before detector acceptance effects) for the 
parameters as used in the generation and at their required values. The ratio of the two 
spectra is used as a weight to modify the simulation shape in equations 7 and 10. To 
estimate the systematic uncertainty, the input value is changed within its error and the 
procedure is repeated. 

A similar approach is employed to allow for the uncertainty of the means (X E (D)) and 
(X E (B)). JETSET is used to generate the Xg distributions of all charm states according 
to (X E (D*)) = 0.510±0.005±0.008, (X b E (B)) = 0.702±0.008 [16]. The energy spectrum 
of D mesons in the B rest frame was measured by CLEO [17]. This spectrum includes 
the contributions from B — > D X and B — > DD X. It can be parameterized in terms of 
a Peterson function with = 0.42 ±0.07 [16]. 

The corrections are applied to all simulated charm ground state hadrons separately 
for bb and cc events. The resulting Xg distribution of the sum of all charm hadron 
ground states in cc events is found to be in agreement with the corresponding average of 
(X C E (D°,D + )) = 0.484 ±0.008 [16]. Here the effect of gluon splitting into cc is taken into 
account. The systematic uncertainties are calculated separately for (X E (D)), (X E (B)) 
and 

The 6-hadron lifetimes are corrected separately for B + , B°, A b and B®. Here the world 
averages t(B°) = 1.56 ± 0.04, t(B+) = 1.65 ± 0.04, t(B° s ) = 1.54 ± 0.07 and r(A 6 ) = 
1.22 ± 0.05 ps [9] are used to correct the simulation. For the systematic uncertainties 
from this source, all the b lifetime distributions are regenerated with a change of one 
standard deviation and the fit is performed again. Similarly the c-hadron lifetimes are 
also corrected separately for D + , D°, A c and D+. Here the values t(D°) = 0.415 ±0.004, 
t(D + ) = 1.057 ± 0.015, t(D+) = 0.467 ± 0.017 and r(A c ) = 0.206 ± 0.012 ps from [9] are 
taken. 

The separation between bb and cc events obtained from the 6-tagging depends on the 
production rates of D + and D° mesons in cc events. The rates of charm hadrons in the 
hemisphere opposite to the reconstructed D are therefore fixed to the present averages 
/(£>+) = 0.221 ± 0.020, f(D+) = 0.112 ± 0.027 and f{c baryon ) = 0.084 ± 0.022 [18]. The 
D° rate is calculated from these numbers according to: 

f(D°) = 1 - /(£>+) - /(£>+) - f(c baryon ) . (17) 

A variation of one standard deviation on each fraction is included in the systematic error, 
leaving the D° fraction free to keep the sum constant. 

The effect due to the efficiency of the 6-tagging was studied in reference [2] using a 
tuning determined independently on data and simulation. A residual difference in the b 
efficiency of 3 % per jet between data and simulation was found. The corrections to the 
physics parameters in the simulation mentioned above account for this difference. This 
systematic error is therefore already included in the physics corrections. Furthermore, 
the effect due to the resolution of the 6-tagging is determined by interchanging the 6-tag 
tunings of data and simulation with each other. 

The rate of gluon splittings into cc pairs is varied by one standard deviation. 

The relative rate of D mesons from bb and cc events is not a free parameter in the 
asymmetry fit. Therefore the ratio is fixed to the DELPHI measurement [19] using the 
same data and varied by one standard deviation. 

The mixing correction for A b FB is discussed in section 6. The systematic error is 
obtained by varying separately each parameter used to obtain the B decay rates into D 



15 



mesons by one standard deviation. The effect of the variation is studied directly on the 
asymmetry fit; this allows for the lifetime dependence of the B° — B° mixing. The total 
error is then calculated taking the correlation between the parameters into account. The 
effect of the oscillation frequency error is small compared to that from the uncertainty of 
the B — > D rates. 

Differences between the signal and background efficiency as a function of cos# are 
considered in the calculation of the probabilities from the simulation. Because the asym- 
metry evaluation depends on the ratio of D to D at a given cos#, the sensitivity to 
efficiency variations (which are largely independent of the nature of the D) is small. The 
systematic error due to the fit method is estimated by comparing the results of the x 2 fit 
to the results obtained assuming a Poisson distribution and neglecting in both cases the 
error on the simulation. The observed difference is included in the systematic error. 

For all decay modes the relative normalization Rs/b is obtained from a fit of the 
simulated D signal and background to the data. A variation of ±5 % (±10 % for the 
off-peak data) is included in the systematic error, not only to account for the error of the 
fitted Rs/b, but also for uncertainties in the agreement of the shape of the mass difference 
signals in data and simulation. 

The contribution of partially reconstructed D decays depends on the efficiency to 
reconstruct such n s i + X combinations, as well as on the total rate of D* + — > D n + 
decays in hadronic Z events. The differences for the background normalizations between 
data and simulation average around 10 % whereas the total rate of D* + — > D°n + decays 
is known at the 5 % level. Combining these, the contribution to the systematic error is 
estimated as a ±15 % variation of the prediction of the simulation. 

The three classes of the combinatorial background (uds, c and b) have a small remain- 
ing asymmetry. The asymmetry of the uds quark background is taken from the simulation 
and subtracted in the fit. The charge correlation for the combinatorial background from 
b and c events is taken from the simulation. The agreement between data and simulation 
is tested using the side bands, where good agreement is found. 30 % of the effect due to 
the background asymmetry is considered in the systematic error on the asymmetry. 

The contributions to the systematic errors for the combined fit of the charm and 
bottom asymmetries are listed in table 8. The relative sign of the systematic error 
indicates the direction in which the results change for a particular error source. 



The results of the 2 parameter fits of the c and b asymmetries for the different energies 
are given in tables 9, 10 and 11. 

Taking the correlations into account, the combination of the results of the different 
samples at the peak energy leads to: 



with a statistical correlation of —0.27. The average centre-of-mass energy is y/s = 91.235 
GeV. In figure 4 the fit results for the forward-backward asymmetries of the different 
samples are compared to the average. The forward-backward asymmetry averaged over 
all samples as a function of cos 9 thrust is shown in figure 5. 



9 Fit results 



A£ B (v/s = 91.235 GeV) 
A b FB (y/s = 91.235 GeV) 



0.0659 ± 0.0094 (stat) 
0.0762 ± 0.0194 (stat) 
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decay mode 



^-fb 



FB 



correlation 


X 2 /N.D.F. 


-0.26 


0.97 


-0.28 


0.87 


-0.28 


0.79 


-0.26 


1.10 


-0.23 


1.26 


-0.27 


0.92 


-0.33 


1.04 


-0.37 


0.80 


-0.33 


0.82 


-0.27 


0.54 



D* + 

D*+ 
D *+ 

D *+ 

D *+ 

D*+ 
D+ 
D° 
D° 



(K 
(K 
(K 
(K 
(K 
(K 
K~ 
K~ 



~7T + )7T + 

^n + n^n + )n + 
"7r + 77)7r + 
~ n + v)it + 
~e + v)it + 

-7T+(7r ))7T + 
7T + 7T + 
7T + 

7T+(7T ) 



0.0590 ± 0. 

0.0738 ± 0. 

0.0255 ± 0. 

0.1130 ± 0. 

0.1201 ± 0. 

0.0720 ± 0. 

0.0567 ± 0. 

0.0431 ± 0. 

0.0534 ± 0. 



0239 
0237 
0298 
0328 
0434 
0209 
0256 
0376 
0405 



0.0498 ± 
0.1125 ± 
0.0914 ± 
0.0131 ± 
0.0218 ± 
0.1037 ± 
-0.0066 ± 
0.0805 ± 
0.0936 ± 



0.0481 
0.0560 
0.0654 
0.0738 
0.0996 
0.0468 
0.1298 
0.0489 
0.0489 



Average 



0.0659 ± 0.0094 



0.0762 ± 0.0194 



Table 9: Results of the two parameter fit to the individual decay modes. The average 
centre-of-mass energy is 91.235 GeV. The x 2 /N.D.F. of the averages is 0.54. 



decay mode 



D* + 

D*+ 

D*+ 

D*+ 
D *+ 

D *+ 

D+ 
D° 
D° 



(K~7Y + 'K~'K + )'K~ 

(i^~7r + 77)7T + 

(K-^+v)ti + 

(K-e + v)ir + 

{K-7l + (7l°))7l + 
K~7T + 7T + 

K--K+ 



Average 



A 



FB 



0.0286 ± 
-0.1893 ± 
-0.2496 ± 
-0.0438 ± 

0.0342 ± 
-0.0514 ± 
-0.0600 ± 

0.2614 ± 

0.1240 ± 



0.1026 
0.0950 
0.1187 
0.1656 
0.1638 
0.0773 
0.0903 
0.1609 
0.1333 



-0.0496 ± 0.0368 



FB 



0.1417 ± 
0.2492 ± 
0.0915 ± 
-0.3322 ± 
0.4078 ± 
0.0475 ± 
0.0158 ± 
-0.0195 ± 
-0.1309 ± 



0.1940 
0.2170 
0.2934 
0.3417 
0.3746 
0.1763 
0.3888 
0.2044 
0.1704 



0.0567 ± 0.0756 



correlation 



-0.26 
-0.27 
-0.31 
-0.38 
-0.23 
-0.27 
-0.26 
-0.34 
-0.34 



-0.28 



Table 10: Results of the two parameter fit to the individual decay modes. The average 
centre-of-mass energy is 89.434 GeV and the x 2 /N.D.F. of the averages is 0.96. 

The results of the different samples at the off-peak energies are combined to give: 



-0.0496 ± 0.0368 (stat) 

0.0567 ± 0.0756 (stat) 

0.1180 ± 0.0318 (stat) 

0.0882 ± 0.0633 (stat) 

The statistical correlation is -0.28 for ^1 = 89.434 GeV and -0.26 for ^1 = 92.990 GeV. 
The results of the fit are shown in figures 6 and 7. The averages of the c and b asymmetries 
for the different energies are shown in figure 8. 



A c FB (^s~ = 89.434 GeV) = 

Afb(Vs = 89.434 GeV) = 

A c fb(Vs = 92.990 GeV) = 

Afb(Vs = 92.990 GeV) = 
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decay mode 



A c 

FB 



A b 

FB 



correlation 



D*+ -> (^-7T+7r-7r+)7r- 

D*+ -> (X-7T + 77)7T + 

D*+ -> (ir-e + z/)7r+ 

Z)*+ -> (^-7T+(7r ))7T + 



0.1090 ± 

0.1281 ± 

0.1778 ± 

0.1190 ± 

0.3113 ± 

0.1049 ± 

-0.0016 ± 

0.2065 ± 

0.0384 ± 



0.0772 
0.0770 
0.1047 
0.1212 
0.1452 
0.0720 
0.0902 
0.1253 
0.1285 



0.1214 ± 
0.2124 ± 
0.1690 ± 
0.2652 ± 
-0.0129 ± 
-0.0798 ± 
0.1387 ± 
0.0049 ± 
0.1227 ± 



0.1616 
0.1765 
0.2535 
0.2655 
0.3003 
0.1494 
0.3652 
0.1575 
0.1543 



-0.25 
-0.27 
-0.29 
-0.32 
-0.28 
-0.25 
-0.19 
-0.33 
-0.34 



Average 



0.1180 ± 0.0318 



0.0882 ± 0.0633 



-0.26 



Table 11: Results of the two parameter fit to the individual decay modes. The average 
centre-of-mass energy is 92.990 GeV and the x 2 /N.D.F. of the averages is 0.50. 



10 The effective electroweak mixing angle 

In order to obtain the effective electroweak mixing angle from the bare asymmetries 
A° FB and A° FB at the nominal Z mass, small corrections have to be applied to the measured 
forward-backward asymmetries. The measured asymmetries for the different centre-of- 
mass energies have to be corrected to ^Js = Mz. The slope of the asymmetry around 
M z depends only on the axial vector coupling and the charge of the final fermions. It is 
thus independent of the pole asymmetry itself. In addition QED corrections, in particular 
initial state photon radiation, reduce the effective centre-of-mass energy. Finally the dia- 
grams from 7 exchange and 7Z interference result in a small correction to the asymmetry. 

The corrections have been determined using the ZFITTER program [20] to be: 



source 


SA e FB 


SA b FB 


v^ = M z 


-0.0034 


-0.0013 


QED effects 


+0.0104 


+0.0041 


7, 7Z 


-0.0008 


-0.0003 


total 


+0.0062 


+0.0025 



where the bare asymmetry is given by A% = A f FB + Ei(SA FB )i. The individual QCD 
corrections for this analysis are included in the quoted measured asymmetries for c and 
b. From the fit results for all energies, the bare asymmetries for c and b are calculated to 
be: 

A% = 0.0715 ± 0.0093 

A°/ B = 0.0793 ± 0.0194 

with a total correlation of —0.22. The quoted error is the combination of the statistical 
and systematic errors at the three different energies. This leads to an effective electroweak 
mixing angle: 

sin 2 6$} = 0.2332 ± 0.0016 . 
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11 Conclusion 

A measurement of the forward-backward asymmetries A C FB and A b FB at LEP 1 energies 
is performed using about 3.5 million hadronic Z decays collected by the DELPHI detector 
in the years 1992 to 1995. The heavy quark is tagged by the exclusive reconstruction of 
D meson decays in the modes D* + -> D°tt + , D+ -> K-tt + tt + and D° -> K'^^). The 
forward-backward asymmetries for c and b quarks at the Z resonance are determined to 
be: 

A c FB (y/s = 91.235 GeV) = 0.0659 ± 0.0094 (stat) ± 0.0035 (syst) 
A b FB (y/s~ = 91.235 GeV) = 0.0762 ± 0.0194 (stat) ± 0.0085 (syst) 

with a total correlation of —0.22. 

The analysis of the off-peak data from 1993 and 1995 leads to: 

A c FB (y/s = 89.434 GeV) = -0.0496 ± 0.0368 (stat) ± 0.0053 (syst) 

Afb(Vs — 89.434 GeV) = 0.0567 ± 0.0756 (stat) ± 0.0117 (syst) 

A c fb(Vs = 92.990 GeV) = 0.1180 ± 0.0318 (stat) ± 0.0062 (syst) 

Afb(Vs — 92.990 GeV) = 0.0882 ± 0.0633 (stat) ± 0.0122 (syst) 

with a total correlation of -0.28 at ^1 = 89.434 GeV and -0.24 at ^ = 92.990 GeV 
respectively. 

The results are in good agreement with other LEP measurements [14,21] using re- 
constructed D mesons. The use of the full available sample of the reprocessed data for 
the years 1992 to 1995 leads to a significant improvement in statistical precision as com- 
pared with the previous results from DELPHI [1]. The results and the obtained energy 
dependence are consistent with the predictions of the Standard Model. 

From the corresponding bare asymmetries for c and b quarks, the effective electroweak 
mixing angle is determined as: 

sin 2 9 l e e $ = 0.2332 ± 0.0016 . 



This result is in good agreement with the determinations of the effective electroweak 
mixing angle from several LEP and SLD measurements [18]. 
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D* + H> (K"7t + )7t + D* + -> (K"7t + (7t°))7t + 




Am [GeV/c 2 ] Am [GeV/c 2 ] 



D* + -> (K"7t + 7t"7t + )7t + D* + -> (K'7t + YY)7t' 




Am [GeV/c 2 ] Am [GeV/c 2 ] 

Figure 1: The mass difference distributions Am for the different decay modes. Am is 
defined as the difference between the mass of the D* + and the D° candidate. The data 
are compared to the simulation. Contributions from reflections, partially reconstructed 
D* + decays (7r s ; + X) and combinatorial background are also shown. See section 5 for 
the discussion of these contributions. 
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Figure 2: The mass difference distributions Am for the semileptonic decay modes are 
shown in the two upper plots. Am is defined as the difference between the mass of the 
D* + and the D° candidate. The data are compared to the simulation. Contributions 
from reflections, partially reconstructed D* + decays and combinatorial background are 
also shown. The lower diagrams show the D + and D° mass distributions. For the D° the 
background distribution for candidates with wrong mass assignments is also shown. In 
the D° case, the peak in the K~ir + mass distribution comes from D° decays into K~tt + , 
and the broad enhancement at lower values is from the K~7t + tt° decay. 
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Figure 3: The X E , tr(V ev ) and Q ■ cos(9 t hrust) distribution for the signal region of the 
D* + — > (K~ir + )ir + decay mode. The data are compared to the simulation where D* + 
from charm and bottom events and combinatorial background are shown separately. See 
section 5.1 for the definition of tr(V ev ). 
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Figure 4: The results of the two parameter fit of the c and b asymmetry at an average 
centre-of-mass energy of 91.235 GeV for the different D samples are shown. The grey 
bands represent the averages over all these measurements. Only statistical errors are 
shown. 
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Figure 5: The c and b forward-backward asymmetries at an average centre-of-mass 
energy of 91.235 GeV as a function of cos 6 thrust- Only statistical errors are shown. The 
bands represent the fit results. 
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Figure 6: The results of the two parameter fit of the c and b asymmetry at an average 
centre-of-mass energy of 89.434 GeV for the different D samples are shown. The grey 
bands represent the averages over all these measurements. Only statistical errors are 
shown. 
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Figure 7: The results of the two parameter fit of the c and b asymmetry at an average 
centre-of-mass energy of 92.99 GeV for the different D samples are shown. The grey 
bands represent the averages over all these measurements. Only statistical errors are 
shown. 
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Figure 8: The c and b forward-backward asymmetries for the different average 
centre-of-mass energies. The SM prediction is also shown (Mz = 91.187 GeV/c 2 , 
m top = 175.6 GeV/c 2 , m^ = 300GeV/c 2 , a s (M|) = 0.120, a(M|) = 1/128.896). 



